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Abstract 



Raman interaction of a trapped ultracold ion with two traveling wave 
lasers has been used extensively in the ion trap experiments. We solve this 
interaction in the absence of the rotating wave approximation by a continued 
fraction, without considering the restriction of the Lamb-Dicke limit and the 
weak excitation regime. Some interesting characteristics of the ion-trap sys- 
tem, particularly for the ion outside the weak excitation regime are found. 
Finally, a comparison of our results with the solution under the rotating wave 
approximation is made. 
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Quantum effects in the motion of trapped ultracold ions have been drawn much attention 
in recent years since the success that a single ultracold ion was confined in the ground state 
of a Paul trap' 1 '. Considerable progress has been made so far, such as the preparation and 
measurement of various nonclassical states' 2 ', and two qubits logic gate operation' 3 ', which 
are of great current interest in ion trap experiment. 

It has been proven that a trapped ultracold ion experiencing laser waves can be described 
by the Jaynes-Cummings model(JCM)' 4 ' under the rotating wave approximation(RWA). 
The JCM, describing a two-level system coupled to a harmonic oscillator under RWA in 
the case of near resonance and weak coupling, has successfully dealt with the cavity QED 
problems, in which a two-level atom interacts with a monochromatic cavity mode with the 
form of harmonic oscillator. In the trapped ion configuration, the ion is simplified to be 
two levels and supposed to move within the region of the space much smaller than the 
effective wavelength of the laser wave(Lamb-Dicke limit), the vibrational level of the trap is 
quantized as the harmonic oscillator, and the laser radiating the ion is generally supposed 
as the classical form of a standing or traveling wave. Different from the situation of cavity 
QED, the strength of the coupling between the ion and the oscillator can be conveniently 
adjusted simply by changing the intensity of the laser. For example, in the strong excitation 
regime' 5 ' 6 ' 8 ' , the coupling constant, denoted by Rabi frequency, is larger than the trapping 
frequency, which makes RWA no longer valid. 

Under RWA, along with the assumption that the ion is confined within the Lamb-Dicke 
limit(LDL) and weak excitation regime(i.e., the Rabi frequency being smaller than the 
trapping frequency), one may easily obtain three simple cases' 4,6 ', that is, carrier excitation, 
red detuning and blue detuning. If excluding the supposition with respect to the LDL and 
weak excitation regime, Vogel and his coworkers gave a standard approach to such a nonlinear 
JC model' 7 '. However, in case the RWA is excluded, the ion trap problem turns to be non- 
integrable, whose solutions have to be obtained completely by the numerical calculation' 6 ' 8 ]. 
Recently, we analytically treated' 9 ' the interaction in the absence of RWA between a standing 
wave laser and a trapped ultracold ion in coherent state representation' 10 ' 11 ' . Although 
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those analytical results are only some particular solutions to the problem, some important 
information neglected in former works due to use of RWA were presented. However, our 
investigation in that work was only under both the weak excitation regime and the LDL. 
Moreover, the solutions in that work strongly depended on the truncation of the series, 
and thereby the deduction and the forms of the analytical expressions went tedious with the 
increase of the number of terms in the series, which makes it hard to obtain the eigenenergies 
and eigenfunctions rapidly at our disposal. Furthermore, it is Raman configurations that 
are usually applied in actual ion trap experiments for the interaction between lasers and 
ions, instead of the simple case described in Ref.[9]. Therefore, in this contribution, we 
choose Raman A-type configuration(shown in Fig.l) for a study, which corresponds to the 
actual process in the NIST experiments' 1 '. To present a more general description for such 
a model, we will exclude the assumption related to both the weak excitation regime and 
LDL by means of some unitary transformations. Moreover, to overcome the difficulty of 
complicated analytical deduction, our solution will resort to a continued fraction, from 
which the eigenenergies of the system and the coefficients of the series of the eigenfunctions 
can be rapidly obtained by numerical calculation. 

Fig.l describes the interaction of a trapped ultracold ion with two off- resonant counter- 
propagating traveling wave lasers with frequencies uj\ and cj 2 respectively. Both laser beams 
are assumed to propagate along the x axis, and the problem is thereby one-dimensional. 
For a sufficiently large detuning 5, the third level \r > may be adiabatically eliminated, and 
we face an effective two-level system, in which the lasers drive the electric-dipole forbidden 
transition \g >^ \e >. In the frame rotating with the effective laser frequency ul (= lo\— 
the dimensionless Hamiltonian of this system is' 1 ' 2 ' 12 ] 

H =^-a z + a + a + ^(a + e ir > £ + (7-6-^) (1) 

where A = (c^o — with u Q being the resonant frequency of the two levels of the 

ion, and v being the trap frequency. Q is the dimensionless Rabi frequency and 77 the 
Lamb-Dicke parameter. Oi (i = ±, z) are Pauli operators, x = a + + a is the dimensionless 
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position operator of the ion with a + and a being operators of creation and annihilation of 
the phonon field, respectively. The notations '+' and '-' in front of ir/x indicate that the 
laser wave propagates towards x and -x directions, respectively, v is generally supposed to 
be much greater than the atomic decay rate(strong confinement limit) for neglecting the 
effect of the atomic decay. 

We first perform some unitary transformations on Eq.(l)' 12 l, and obtain 

H 1 = UHU + = ^-a z + a + a + g(a + + a)(a + + cr_) + e(o+ + <7_) + g 2 (2) 
( D+{(5) D{(5)\ 

where U = -4= e ina+a / 2 with D((5) = e ^ {a++a ^ 2 , g = rj/2 and e = -A/2. 

\-D+(f3) D(P)J 

Comparing with Ref.[9], Eq.(2) can be regarded mathematically as the simplified form of 
Eq. (6) in Ref . [9] without the two-phonon processes, whereas the optical resonance frequency 
is replaced by the Rabi frequency. Similar to Ref. [9], we set a + — > a, a — > and rewrite 
Eq.(2) in coherent state representation, that is, 

H = % z + a^- + g(a + -^)(<r+ + <r_) + e(a+ + cr_) + # 2 (3) 

where a is a complex number and / exp(— |o;| 2 )|a* >< a*\ = 1. The eigenfunction 
of the Schrodinger equation of Eq.(3) can be supposed to be ^(a) = I I , and the 



eigenenergy is E. Thus we have 



# 2 (a) 



(a + 9)^-^1 = (E - g 2 - ga- e)®! - ^$ 2 , (4) 
da 2 

(a-5)^$2 = (£-£ 2 +0a + e)$2-^$i (5) 

where we have made a transformation <3>i = ^i(a) + ^(a) and <3> 2 = ^i(a) — ^(a)- To 
solve Eqs.(4) and (5), we order £ = a + g, and <3>i = exp (— Combining Eq.(4) with 
Eq.(5) will yield a second-order differential equation 

z(t - 2 ^^ + t 2 ^ - 2 ^ - 2 ^ e + 1 1 + 4 ^ 2 - - 2 ^ 2] ^ 0+ 

[£ 2 - e 2 - ^ - Ag 2 E + 4<? 2 e + 2g(E - e)£]0 = 0. (6) 



By supposing to be a series form = I]^L C n £ n , Eq.(6) becomes a recurrence relation 
about the coefficient C„ 
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a n C n+ i + (3 n C n + 7„C n _i = (7) 

where a n = 2g(n + l)(E-e-l-n), (3 n = [n 2 + 2n(2g 2 - E) + E 2 -e 2 - - 4g 2 E + Ag 2 e] and 
7„ = 2g{E— e— n+1 ). Obviously, by cutting the series at C_i = 0, we have C\/Cq = — 0o/ a o, 
and other coefficients C n (n = 2, 3, ■ ■ •) will be denoted by C . Making a slight transformation 
for Eq.(7), we obtain following continued fraction, 



' n ~~f> 7n-l Q n-2 ~f> Qra+17n+2 V / 

Pn-1 ~5 ' 72"! Hn+l T, "n + 2')'n + 3" 
Pn~2" 5l _ 7iao// 3 /^n+2 



from which the eigenenergies E can be easily obtained, and the coefficients C n (n — 1, 2, • • •) 
can also be solved' 13 ]. 

Prior to the numerical calculation for Eq.(8), we first consider a special case of the 
problem, that is, the case of LDL and large detuning in the strong excitation regime. At 
this time, Eq.(3) reduces to the mathematical form similar to the carrier excitation, i.e., 
H = Qa z + a-^ + e(er + + cr_). With above procedure, we can obtain 
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p n = 0^ E = n± ] J— + e 2 

which is in good agreement with the solution in Fock state representation. Except for this 
special case, the other situations of the problem have to be solved numerically from Eq.(8). 
Obviously, the convergence of the coefficients C n is the necessary condition in the calcu- 
lation. As this condition has been discussed in Ref.[13], we will not repeat it here. The 
numerical calculation of Eq.(8) demonstrated in Figs. 2, 3 and 4 may present us following 
interesting results: 

(i) the eigenenergies of the system increase with the enhancement of Q and 77, which means 
that it is more difficult to stably confine the ultracold ion outside the LDL and weak exci- 
tation regime; 

(ii) the increase of the eigenenergies with the enhancement of r] is suppressed by the large 



detunings. In view of physics, the only reasonable explanation for this case is that the in- 
teraction between the laser and the ion plays main role in the system. With the increase 
of the detuning, this interaction goes weaker and weaker. Therefore, although the energy 
of the ion itself is enhanced with the increase of 77, the energy of the total system remains 
nearly constant; 

(hi) Fig. 4 presents the relation between Q and A for certain values of E and rj. From the 
forms of the curves, we may guess that the relation among E, A, 77 and Q is of the form of 
= f(E)[l — h(r/)A 2 ] with f(E) and h(rj) being positive definition functions of E and 77, 
respectively. 

Moreover, from the specific calculation, we found that, for a certain Q, with the increase 
of rj, the convergence of C n would get worse and worse, particularly for the large detuning 
case. But the larger the value of Q, the larger the value of rj for satisfying the convergence 
of C n . The problem with respect to the convergence will be discussed physically later. 

Returning Eqs.(4) and (5), we can find that $1 = exp (g£)<f> is another solution for 
£ — a + g. Repeating above deduction, we will obtain Eq.(8), whereas in this case, a n = 
2g(n+l)(E-l-e-n), (3 n = (n 2 -2nE -4ng 2 + E 2 -e 2 - -4g 2 ) and 7ri = 2g(n-e-E). 
Similar to above procedure, we studied the relation among E, A, Q and 77, as shown in 
Figs. 5, 6 and 7. Comparing with Figs. 2, 3 and 4 respectively, we may find that the values 
in Figs. 5, 6 and 7 are slightly different, and their difference become more and more obvious 
with the increase of 77 and f2, or E and 77. Particularly for the case of f2 = 2 and A = 2.0 in 
Fig. 6, the curve no longer upgrades with the increase of 77. 

As we obtained the solutions of $1, $ 2 can be easily solved from Eqs.(4) and (5) with 
the form of $ 2 = En=o( n - E + e)C n (a + g^e^+a). Therefore, the eigenfunction of 
the system is 

1 + .(T£#C n [l-l(n-E + e)](a + gr\ 
*(<*)= = 7:e ±9{a+g) " (9) 

V Ma) J 2 V C n [l + £ (n - E + e)] (a + g)» ) 

From Refs. [9,10,15], we know that the coherent states can be conveniently transformed to 

the Fock states with the formula 
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n >= / dada exp(-|a| 2 )|a* > —La". (10) 
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By means of Baker-Campbell-Hausdorff formula, we may find e ga a n in coherent state rep- 
resentation correspond to the displaced Fock state \n,g >. Therefore, ty(a) in Eq.(9) is in 
fact a superposition of several displaced Fock states. However, to obtain the eigenfunctions 
in the original representation, we should consider the transformation made in Eq.(2). As an 
example, we may treat the special case H = -^a z + a-^ + e(<7+ + er_), whose eigenfunction is 

* s (*) = )=U E<?»« B (ii) 



where A = 1 ± + e 2 + ff , 5 = 1 =f + e2 + If- Then the eigenfunction in the 

original representation is 



1 + — 



\n > 

(12) 



e 2 



where we have made a n to be Fock state \n > by means of Eq.(10), and simply chosen the 

coefficients related to a n to be 1 and a m (m ^ n) to be zero. From above equation we know 

that, in the weak excitation regime Q <^ e, or in the carrier excitation, i.e., e = 0, ^ so (a) oc 
1 \ 

\ \n >, which is identical to the standard solution from Fock state representation^. 
±1/ 

In general, in the strong excitation regime with large detuning, the probability amplitudes 
of up and down states will be different, as shown in Eq.(12). 

As referred to previously, the problem of eigenfunction and eigenenergy for a trapped and 
radiated ion was treated numerically in Ref.[6], in which the Hamiltonian was diagonalized 
in the space spanned by the dressed Fock states. As the solutions of the non-RWA JC 
model are with the forms of coherent states or squeezed states, etc' 9,14,15 ', the accurate 
eigenenergies and eigenfunctions in Ref. [6] might not be in principle solved unless the space 
for diagonalization was spanned by the infinite dressed Fock states, particularly in the case 
of the large detuning and strong excitation regime. That is the reason why Ref. [6] only 
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presented the solutions in the weak excitation regime for the case of near-resonance. In our 
treatment, however, the numerical solution is made based on the continued fraction, which 
is obtained in the coherent state representation. So our solution is still valid even for the 
large detuning and strong excitation regime. 

Before ending our discussion, it is interesting to make a comparison between our solution 
with those under RWA. In general, we may solve Eq.(l) by means of nonlinear JC model as 
proposed in [7] for the ion governed outside the weak excitation regime. However, to make 
the comparison easily and obviously, we start our comparison procedure from Eq.(2). By 
performing a unitary transformation exp[— i(^a z + g ^f L )t\ on Eq.(2), we can obtain E„ = 
+ \if ± \yj4rj 2 (n+ 1) + 1 (n = 0, 1, 2, • • •) in the Fock state representation for ft = 2 
under RWA(see Figs. 8). From Figs. 3, 6 and 8, we see that, for small values of rj, the case 
of n = 2 under RWA is similar to those of Q = 2 in Figs. 3 and 6, where E+ and E~ are 
similar respectively to the very large and zero detuning cases without RWA. Although both 
the solutions with RWA and without RWA include multi-curve, the difference is obvious: as 
it is merely a particular solution, the solution without RWA approaches to the solution of a 
certain excited state under RWA. On the other hand, as RWA merely retains some resonance 
terms in Hamiltonian, the solution under RWA is irrelative to the variation of A and just 
related to a specific value of Q. Moreover, as they are obtained from the solution of Eq.(8), 
the results without RWA are restricted by the convergence of the coefficients of the series. In 
contrast, there is no restriction mathematically for the expressions of the solution under RWA 
even in the case of rj — > oo. However, in physics, rj is impossible to approach to infinity for a 
ultracold ion system, and it must be cut off at a certain value. Therefore, in this sense, the 
convergence condition in our solution of the continued fraction provides a physical restriction 
on the final results, which plays the same role as the restricted condition for the solutions in 
Ref. [9]. Furthermore, further deduction with the transformations exp[— i(^a z + g ^ L )t] and 
exp[— i(^a z + ^)t] performed on Eq.(2) respectively can tell us that, with the increase of 
Q and rj, the difference between the solutions with RWA and without RWA will be larger 
and larger. 
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In summary, we have presented a continued fraction solution to a trapped and laser- 
radiated ultracold ion experiencing two traveling wave lasers in a Raman process in the 
absence of RWA. Instead of the complicated and tedious analytical forms in Ref.[9], we can 
solve the problem rapidly from the continued fraction and present more complete information 
of the system, such as the relation among E, Q, A and rj. As limitation of weak excitation 
regime and LDL as well as RWA is excluded in our treatment, our solution is not only exact, 
but more general. It can be used to investigate quantum properties of the ion-trap system 
in a wider range of parameters, and serve as a comparison with other approximate works, 
as we did in Ref.[16]. However, as Eq.(2) represents a kind of typical non-integrable system, 
our solutions are still only some particular ones. While in comparison with former solutions 
under RWA, we still found some interesting results, particularly for the case outside the weak 
excitation regime and LDL. We should emphasize that, it is hard to carry out the trapped 
ion-laser interaction experimentally outside the weak excitation regime at present due to 
the requirement of very high intensity of the laser, and the laser cooling to the ground state 
of the trapping potential has not been reported yet in this regime. Therefore, we hope the 
solution in our work would be helpful for the future exploration of the quantum properties 
of the ion-trap system in this respect. 

The work is partly supported by National Natural Science Foundation of China under 
Grant No. 19904013. 

Note added: after accomplishing this paper, the author becomes aware of that a workt 17 ! 
related to the fast quantum gate for cold trapped ions has been carried out based on Ref. [12] 
and the treatment similar to our RWA treatment. 
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Captions of the figures 

Fig.l Level scheme of the internal structure of the trapped ultracold ion, where \g ><-> 
|e > is dipole forbidden. 

Fig. 2 Variation of E with respect to A in the case of $i = exp (— g^)<p, where the 
dotted, dashed, solid and dash-dotted curves represent rj =0.2, 0.4, 0.6 and 0.8, respectively. 
For clarity of the figure, we only plot the cases for Q =2, 4 and 6. 

Fig. 3 Variation of E with respect to rj in the case of <3>i = exp (— g£)0, where the solid, 
dotted, dash-dotted and dashed curves represent A =0.2, 1.6, 2.0 and 3.0, respectively. For 
clarity of the figure, we only plot the cases for Q =2, 4 and 6. 

Fig. 4 Variation of Q with respect to A in the case of <3>i = exp (— g£)0> where the solid, 
dashed and dotted curves represent rj =0.02, 0.4 and 0.6, respectively. We only plot the 
cases for E =3 and 5. 

Fig. 5 Variation of E with respect to A in the case of $i = exp (g£)0- The rest are the 
same as that in Fig.2. 

Fig. 6 Variation of E with respect to i] in the case of $i = exp (g£)0- The rest are the 
same as that in Fig. 3. 

Fig. 7 Variation of f2 with respect to A in the case of $i = exp (gi)4>. The rest are the 
same as that in Fig.4. 

Fig. 8 Variation of E^ with respect to r\ for Q = 2 for the solution under RWA. See 
text. 
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